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We present a search for ten baryon-number violating decay modes of A hyperons using the CLAS 
detector at Jefferson Laboratory. Nine of these decay modes result in a single meson and single 
lepton in the final state (A —>■ ml) and conserve either the sum or the difference of baryon and 
lepton number (BEL). The tenth decay mode (A —>■ p7r+) represents a difference in baryon number 
of two units and no difference in lepton number. We observe no signihcant signal and set upper limits 
on the branching fractions of these reactions in the range (4 — 200) x 10“^ at the 90% confidence 
level. 

PACS numbers: ll.30.Fs, 13.30.Ce, 13.30.Eg, 14.80.Sv, 25.20.Lj 


I. INTRODUCTION 

The Standard Model (SM) of particle physics [IH3] 
has had great success in interpreting and predicting ex¬ 
perimental results since its conception in the late 1960s. 
There are, however, features of our universe that are in¬ 
consistent with the SM framework. Astronomical obser¬ 
vations suggest that our universe is dominated by mat¬ 
ter over antimatter n [5]. Sakharov proposed in 1967 
that this asymmetry suggests fundamental interactions 
that violate CP-symmetry and baryon-number conserva¬ 
tion [5]. The observed quark-sector CP violation, com¬ 
bined with baryon-number violating (BNV) processes 
that are allowed by the Standard Model, are insuffi¬ 
cient [7] to account for the observed matter-antimatter 
asymmetry in our universe. A possible explanation for 
this discrepancy is that there are yet-unobserved inter¬ 
actions that violate baryon-number conservation. 

Baryon-number violating reactions are features of sev¬ 
eral theoretical extensions to the Standard Model, per¬ 
haps most notably the SU{5) Grand Unihed Theory 
(GUT) of Georgi and Glashow [8], SU{5) being the larger 
gauge group in which the Standard Model’s SU{3) x 
SU{2) are embedded [9]. The SU{5) theory proposes 
the existence of two new gauge bosons, the X and Y lep- 
toquarks, so called because they allow vertices such as 
q -T XI, where q is a quark and I a lepton. Other ex¬ 
periments have been performed to search for BNV pro¬ 
cesses in decays of the nucleon [IDHIS, r leptons [Sl- 
[I5], top quarks |16j . hadrons with bottom and charm 
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quarks [miiH], and the Z boson m, but no signal has 
yet been observed. The most stringent limits on such 
processes come from nucleon decays, and these have been 
used to constrain BNV decays in higher-generation {i.e., 
c, b, and t) quarks [20]. However, multiple amplitudes 
can contribute to a given decay process and these am¬ 
plitudes can interfere either constructively or destruc¬ 
tively, depending on their relative phases. The theo¬ 
retical calculations that constrain these BNV processes 
do not take into account interference between the ampli¬ 
tudes due to the large parameter space. This allows for 
non-observation in one mode {e.g., decays involving the 
u or d quark) while still being consistent with observ¬ 
able BNV processes in some other mode {e.g., coupling 
to another quark flavor). 

Here we present a search for baryon-number (B) and 
lepton-number (L) non-conserving decays of the A hy- 
peron as a direct probe of couplings of BNV interactions 
to the strange quark. By summing the branching frac¬ 
tions and experimental uncertainties (in quadrature) of 
the six observed A decay modes [12], we find the total 
branching fraction to be 1.001 ± 0.007, implying that 
there is room for yet-unobserved decay modes. 

We investigated eight decay modes in which the A de¬ 
cays to a charged meson and a charged lepton, conserving 
charge in all decays. The meson is either a or and 
the lepton is either a or a We produced the A 
by means of a photon beam incident on a liquid hydro¬ 
gen target through the exclusive reaction ■jp —>■ K'^A. 
For these eight modes we can completely reconstruct the 
three final state particles. We also searched for the de¬ 
cay of a A to a Kg and a neutrino, which must be in¬ 
ferred from the missing momentum. Selection of these 
nine channels is motivated by searching for decay of the 
A to a lighter pseudo-scalar meson. In each case the final- 
state meson is included to preserve charge and angular 
momentum conservation. Thus, violation of L in these 
reactions is a consequence of B violation rather than a 
primary motivation. 

In addition, we searched for the BNV decay of the A to 
an anti-proton and a 7r+, for which we can completely re- 
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Decay 

AB 

XL 

A{B-L) 

detected 

A ^ K+£- 

-1 

-fl 

-2 

K+,K2,£- 

A K-£+ 

-1 

-1 

0 

K+,K-,£+ 

A —>■ 

-1 

-fl 

-2 

K+,ti+,£- 

A —>■ 'k~£'^ 

-1 

-1 

0 

K+,n-,£+ 

A ^ K^u, K%v 

-1 

±1 

0, -2 

, TT^, n~ 

A —>■ 

-2 

0 

-2 

K+,p, Tr+ 


TABLE I. Properties of conserved quantities for each A BNV 
channel. Here I represents a lepton, either e or /r. The right¬ 
most column shows the detected final-state particles for each 
channel. Subscripts denote recoil (r) or decay (d) kaons. 

construct the final state. The A —p7r+ reaction presents 
an opportunity to search for A-A oscillations, i.e. a pro¬ 
cess by which the A oscillates into its anti-particle coun¬ 
terpart (A), which then undergoes a standard-model de¬ 
cay (A —>■ p7r+). Though there is not a simple way to pic¬ 
ture this reaction proceeding via X boson coupling, theo¬ 
retical and experimental {e.g. [211122]) work has been per¬ 
formed by other groups looking for similar oscillations of 
the neutron. Such baryon-antibaryon oscillations would 
have far-reaching implications and are often held up as 
evidence for high-energy theories ranging from see-saw 
models (231 [23] to extra dimensions [2S]. 

The properties of these reactions are summarized in 
Table |T| These specihc decays were chosen for several rea¬ 
sons: 

• Each reaction shows evidence of AS yf 0 and/or 
XL = ±1. 

• Each reaction conserves electric charge and angular 
momentum. 

• This selection includes reactions that either pre¬ 
serve or violate B — L, a, conserved quantity pro¬ 
posed by several GUTs [26] . 

• The CLAS detector is optimized to reconstruct the 
final-state charged particles produced in each re¬ 
action, except for the neutrino which can be in¬ 
ferred by calculating the missing 4-momentum in 
the event. 


Analysis overview 

This analysis was performed in three stages. Here, we 
present a brief outline; details of each stage will be given 
in the following sections. 

A identification. In order to assess the sensitivity 
of our study, we first determined the number of A hy- 
perons produced during the run period. We did so by 
considering the charged decay mode, A — pTr~. We ap¬ 
plied a set of simple cuts on kinematic observables and 
timing for the recoil K~^ to effectively identify potential 
A events. We used this sample to determine the total 
number of "fp —>■ K'^A —>■ K'^ptt~ events detected, and 
then acceptance corrected to find the total number of 


7 P —>■ K'^A —> K'^pir events produced during the run 
period. 

Channel-specific tuning. When searching for evi¬ 
dence of the BNV decays listed in Table [l] we performed 
a blind search. We developed a set of background sepa¬ 
ration cuts for each BNV channel, based on timing infor¬ 
mation for all charged final-state particles and kinematic 
observables for the event. In developing these cuts, we 
balanced the optimization of kinematic cuts for both dis¬ 
covery and upper-limit sensitivity by maximizing a fig¬ 
ure of merit (approximately the BNV signal efficiency 
divided by the square root of the number of background 
events). We assessed the signal efficiency using a Monte 
Carlo technique and the background size using side-bands 
of the blinded signal region. This step provided a set of 
cuts for each BNV channel that reduces background and 
provides optimal analysis power. 

Unblinding. We then unblinded the signal regions 
of kinematic variable plots, and determined whether 
a signal is present. For nine of the decays, the ex¬ 
pected backgrounds are 0 or 1 event and so we used the 
Feldman-Cousins method m to determine upper bounds 
on branching fractions; for the remaining channel, the 
backgrounds are higher and so we scanned the relevant 
parameters in our fit to determine the 90% coverage for 
the number of signal events in the dataset. 


II. THE CLAS DETECTOR AND DATASET 

The CLAS detector is described in detail else¬ 
where [2H]- The dataset comes from the gll run pe¬ 
riod, which collected data during May and June, 2004. 
A bremsstrahlung photon beam was produced by a 4.023- 
GeV electron beam incident on a gold radiator. Electrons 
were provided by CEBAF (Continuous Electron Beam 
Accelerator Facility) in 2-ns bunches. Photon energy and 
timing information were provided by a tagging spectrom¬ 
eter which directs the electrons after the radiator through 
a magnetic field and onto a set of scintillators, providing 
photon energy resolution of 4.0 MeV. 

For this analysis we made use of the CLAS drift cham¬ 
ber and toroidal magnet systems to measure the mo¬ 
menta of the charged final-state particles. Velocity mea¬ 
surements are made by a start counter (consisting of 
scintillators placed within 11.6 cm of the target) and a 
set of time-of-flight (TOF) scintillators located approxi¬ 
mately 4 m from the target. Timing information from the 
photon tagger is combined with this system to calculate 
a velocity for each charged track and this is compared 
with the particle hypothesis in the particle identification 
(PID) algorithm. More details of the gll dataset, cali¬ 
bration procedures, and systematic studies can be found 
in {e.g.) [22]. 
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III. IDENTIFICATION OF ~fp K+A K+pn- 
EVENTS 


tern, and c is the speed of light. This PID cut (and later 
cuts) is based on the difference of these two quantities, 


In order to compare the branching fraction for a BNV 
decay mode to that of the standard-model A —)• pTT~ 
decay, we must first assess the number of standard-model 
decays that occurred during the data-taking period. We 
did so by investigating the (exclusive) jp —>■ A+A —>• 
K~^pTr~ reaction. Earlier studies [29] have shown that the 
7 P —>■ AT+A signal is easily separable from background 
when all three final-state particles are reconstructed. 

Because a different set of background separation cuts 
will be applied to the BNV channels, we must correct 
the number of reconstructed signal events (Vj-ec) to find 
A^prod, the number of yp —>■ AT+A —AT+ptt” events that 
were produced in CLAS during the run period. In order 
to calculate the efficiency of the detector and analysis 
cuts, we generated 3 x 10® Monte Carlo (MC) yp —>■ A'+A 
events, and weighted the distribution of these events in 
cos9^^ (the Ar+ production angle in the center-of-mass 
frame) according to published da/d cos 9 data [53]. Pho¬ 
ton energies between 0.909 GeV (threshold for A+A pro¬ 
duction) and 3.860 GeV (the upper-limit of the pho¬ 
ton tagger during data taking) were generated from a 
bremsstrahlung spectrum given incident electrons of en¬ 
ergy 4.023 GeV (matching the gll run conditions). We 
then used the collaboration-standard GEANT-based [50] 
software suite (GSIM) to model the GLAS acceptance, 
allowing GEANT to produce the A —>• p7r“ decay ac¬ 
cording only to phase space constraints. 

We began the data reduction by selecting from the 
dataset all events in which three reconstructed final-state 
tracks, two of positive charge and one of negative charge, 
were coincident with a tagged photon. To each event, we 
assign the mass hypothesis consistent with A'”*', p, and 7r“ 
final-state particles, selecting the permutation of positive 
tracks as that with the value of the invariant mass of the 
p and 7r“ candidate tracks, INV{p,tt~), closest to the 
nominal A mass. We then apply a loose cut on the square 
of the total missing mass, of each event, keeping 

events for which MM'^ G [—0.002,0.0005] GeV^/c'^. 

To the hypothetical Ar+ track we applied a PID cut 
based on timing and momentum information from GLAS. 
GLAS measures the time-of-flight, tofm, he. the time 
elapsed between the primary event vertex and the track’s 
triggering of the TOE scintillators. We also calculate 
a hypothetical time-of-flight, tofh, based on our mass 
hypothesis and tracking information: 


tofh = 


_d 

/3c 


d 

c 


1-k 




( 1 ) 


where p is the track’s momentum (determined by its cur¬ 
vature through the magnetic field of GLAS); /3 is related 
to the track’s velocity (determined from tracking and 
timing information), j3 = v/c] ruh is the hypothetical 
mass for the track (493.7 MeV/c^ for the Ar+); d is the 
path length of the track from the vertex to the TOE sys- 


Atof = tofh - tofm- (2) 

In the case where the mass hypothesis is correct, we ex¬ 
pect the measured and hypothetical tof to be roughly 
identical (modulo timing resolution), and thus Atof « 0. 
If the hypothetical mass is greater (less) than the parti¬ 
cle’s actual mass, then we expect the Atof to be greater 
(less) than zero. We found that suitable separation of 
Ar+ candidate tracks from non-Ar+ tracks is achieved by 
a two-dimensional cut, keeping events for which 

|Ato/| < (1.8 ns) exp ^ -k0.15 ns. (3) 

The Atof versus |p| plane for signal Monte Garlo events 
is shown in Eig. [^ 

We then make a further cut on the missing mass off of 
the Ar+, 

MM{K~^) = {Pj+Pt-PK^, (4) 

where p^^ pt, and px are the four-momenta of the in¬ 
cident photon, target proton, and recoil respec¬ 

tively. We kept events for which MM(A'+) is in the range 
[1.05 GeV, 1.165 GeVj. In addition to non-strange back¬ 
grounds, this cut removes contamination from photopro¬ 
duction of higher-mass hyperons that include A in their 
decay chain (predominantly yp —>■ AT+E® —AT+yA). We 
also applied geometrical fiducial cuts, omitting events 
from regions of the detector for which our simulation is 
inaccurate. After all cuts, we found that data and MG 
distributions are quite similar, as shown in Fig. [^ 

After these cuts, we identified the signal events by 
inspecting a histogram of INV{p,tt~). Fig. shows 
that this distribution is exceptionally clean. We modeled 
these data with a third-order-polynomial for background 
processes and a double-Gaussian for signal processes. We 
then extracted the number of signal events by taking the 
excess of the data histogram above the background func¬ 
tion in all of the histogram bins within 0.015 GeV/c^ of 
the nominal A mass, yielding Vj-ec = 1.861 x 10® recon¬ 
structed A —p7r“ signal events. Because the shape and 
magnitude of the background and the magnitude of the 
signal are dependent on kinematics, we vetted the above 
estimate by separating the data into ten bins in cos 9^m 
and performing the fit in each bin. This method again 
yields an estimate of = 1-861 x 10®. 


Acceptance correction 

With the number of reconstructed yp —>■ AT ■''A —>■ 
K'^pTT~ events in hand, we could then correct for the 
effects of the detector’s acceptance and the efficiency of 
our analysis cuts to estimate A^prod- We began by ap¬ 
plying the cuts described above to the standard-model 






IV. BNV CANDIDATE SELECTION AND 
OPTIMIZATION 
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Because of the sensitive nature of the possibility of 
BNV discovery, we pursued a blind analysis. For each of 
the BNV A —)■ AB decay modes under investigation (see 
Table |T]), we tuned a set of cuts, C, using the kinematic 
quantities MM(K^), INV{A,B), and MM^, as well as 
the timing information for A and B. As before with 
the Standard Model A decay, we identified BNV signal 
using the INV{A, B) spectrum (except for the A — K^v 
channel described below). In tuning each set of cuts, 
we tried to strike a balance between reducing the large 
number of non-signal events and maintaining acceptance 
for a potentially small BNV signal. 

Punzi |31| has proposed a figure of merit for performing 
such optimizations, which has been used in several other 
searches for rare reactions. For a set of cuts C the Punzi 
figure of merit, 'P{C), is defined as 

no - - 

62 2a/B(^ -k 6^62 + AB{C) 

(7) 

where e{C) is the efficiency of the cuts when applied to 
signal, B{C) is the number of background events pass¬ 
ing cuts C, and a and b are the number of standard 
deviations corresponding to the analysis-defined signif¬ 
icance and statistical power. In this analysis, we chose 
6 = a = 4, indicating a 4cr confidence level. With these 
choices, the figure of merit simplifies to 


FIG. 1. (Color Online) Atof versus momentum (p) for 
candidate tracks for data (top) and MC (bottom) that pass 
the cut on MM’^. The red curves show the upper and 
lower bounds of the A"*" PID cuts. Several additional struc¬ 
tures are apparent in data distribution. The band of events 
with Atof ^ 1 ns represents from decay, that with 
Atof < 1 ns represents protons, and the less-populated bands 
represent events matched with photons from a different 2-ns 
beam bunch. The majority of these are removed by further 
cuts. 


Monte Carlo (MC). Because the generation of MC events 
matches the data in photon energy distribution and kine¬ 
matics, we separated the data and MC coarsely into ten 
bins in cos0^. In each bin, we calculated the accep¬ 
tance CpT^- (cos t>y simply dividing the number of MC 
events that pass CLAS simulation and acceptance cuts 
by the total number of MC events generated. We used 
this factor to determine the number of A+A —>■ A+p7r“ 
events produced during the run in each cos 9^^ bin, and 
completed the calculation by summing over the angular 
bins: 


A^prod — ^ ^ 

COS 


A^rec(cOs6lg„) 

ep^-{cos 


= 3.71 X 10^ 


(5) 

( 6 ) 


V{C) 


e{C) 

16 (2 + ^/B(^) ' 


( 8 ) 


We have tuned our cuts by simultaneously assessing 
e{C) using MC BNV events and B{C) from side-bands of 
the blinded signal region. In all plots that would identify 
any BNV signal {e.g., INV{A,B)) we blind the signal 
region. We postpone the unblinding of the signal regions 
of all of our data plots until after the optimization of 
the analysis cuts, once we are confident that all cuts and 
systematic effects are understood. 

For each BNV A —>■ AB decay under investigation, 
we generate 10® jp —)■ A+A —)■ AB Monte Carlo 
events, matching photon energies to the run conditions 
and A+A kinematics to the measured da/dcos9 (as for 
the standard model MC). We generate kinematics for 
the A —>• AB decays according only to phase-space con¬ 
straints. In the case of the A —)■ A®jz reaction, the sub¬ 
sequent A® —>• TTTT decay is modeled by GEANT at the 
time of detector simulation. 

For each channel, we choose the positive track mass 
hypothesis that yields a value of INV{A,B) nearest to 
the nominal A mass (as was done for the standard-model 
analysis above). The analysis cuts for each channel begin 
with a PID cut on Atof (see eq. for each final-state 
particle. For each particle type, we apply a loose two- 
sided cut in the Atof vs p plane, similar to that applied 
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FIG. 2. (Color Online) INVip,^-) (left), MM{K+) (center), and (right) distributions for standard-model data (filled) 

and MC events (unfilled). In the INV{p,n~) histograms the distributions are of events that pass all cuts, and the vertical 
lines show the limits used in signal counting. In the MM{K^) histograms, the distributions are of events that pass all but the 
cut on MM{K'^) (thus, the peak is still visible in data); the vertical lines show the limits of the cut on MM{K'^). In the 
MM^ histograms, the distributions are of events that pass all cuts. 


to the track for the SM decay. CLAS resolution 
allows us to make these cuts loose; however, the charac¬ 
teristic decay length of is on the same order as the 
dimensions of CLAS and the non-trivial fraction of 
that decay in the detector results in these cuts reducing 
the data sample by approximately half (depending on the 
number of charged kaons detected for each channel). 

A. Example: A — 

Here, we demonstrate this process with the A —> fj.~ 
channel; other channels, with the exception of the K^v 
decay, are analyzed using the same observables. 

For the three charged tracks in each event, we must 
first decide which positive tracks correspond to the recoil 
Ar+ and 7r+. We make this assignment by calculating 
the invariant mass of each positive-negative track pair, 
and choosing the assignment that gives a value nearest 
the nominal A mass. We then apply PID cuts based 
on the Atof method to each of the three particles, the 
boundaries of which are shown in Fig.[^ These cuts are 
loose, and we have found that their efficiencies for each 
particle type are similar when applied to broader MC and 
data samples. 

After the PID cuts have been made, we turn to analysis 
cuts based on kinematic observables. For the A —7r+/j,“ 
channel, we make a symmetric cut on (centered 

about 0 with width wi) and MM{K^) (centered about 
1.1186 GeV/c^, with width W 2 )- In order to find the 
widths, wi and W 2 , which optimize the Punzi metric, we 


uniformly sample forty values for each width, resulting 
in 1600 distinct pairs with 

G [0,0.001] GeVVc'^ (9) 

W2 G [0,0.03] GeV/c^ (10) 

We apply these cuts to both signal MC and data, and 
inspect the resulting , p.~) histograms with the 

signal region of the data histogram blinded. We de¬ 
fine the signal region to be values of INV{'K^,pr) 
within 0.03 GeV/c^ of the nominal A mass (ie., 
[1.086,1.146] GeV/c^), and the side-band regions to be 
within 0.15 GeV/c^ of the peak (excluding the signal re¬ 
gion) . 

We then use the signal MC distributions before and af¬ 
ter cuts to determine the signal efficiency, e{C). We apply 
a simple side-band technique to the data INV{'K^,pr) 
histogram to extrapolate B{C), the expected number 
of background events in the blinded signal region. We 
then use these values to calculate V for each width pair 
(see Fig. [^, and select the width pair that maximizes 
V as optimal. For the A ^ p,~ channel, we find the 

optimal widths to be Wi = 3.25 x 10“^ GeV^/c^ and 
W 2 = 9.00 X 10“^ GeV/c^. The plots in Fig. [^illustrate 
the resulting cuts and the blinded signal plot for this 
channel. 

Fig. [5] demonstrates that the optimal cut scheme for 
this channel is quite restrictive; only five data events pop¬ 
ulate the (blinded) INV{tt~^, p,~) histogram, and none of 
these falls within the side-band regions. We thus estimate 
the expected number of background events in the signal 
region to be zero. By studying the signal MG for this 
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FIG. 3. (Color Online) Fit to the INV{p,Tr~) distributions 
for all events passing the standard model decay identification 
cuts. The lower plot shows the same distribution as the upper, 
but with limited vertical scale. The red curve shows the full 
(signal and background) fit function; the blue curve shows 
only the polynomial background. The magenta lines show 
the limits of integration for the event count. 


channel, we estimate that the efficiency of these cuts, 
e{C), is 7.91% (including the effects of detector accep¬ 
tance). The effects of detector acceptance and analysis 
cuts are similar for all of the charged decay channels. 
Table [IT] lists the properties of the optimal cuts for each 
channel. 


B. A —>■ PTT^ 

Our analysis of the A —>■ p7r+ channel proceeds nearly 
identically to that of the charged A mi channels. We 
first assign the positively charged track mass hypothe¬ 
sis by comparing INV{p,Tr'^) for the two possible track 
permutations. We then apply a PID cut in the Atof vs 
momentum plane for each particle. Threshold p photo¬ 
production via the process yp —?► ppp occurs at a pho¬ 
ton energy of 3.751 GeV; the maximum tagged photon 
energy in our dataset is 3.86 GeV. Unlike in the other 
charged decay channels where there are signihcant num¬ 
bers of each particle type present in the data, the null 
hypothesis suggests that there should be relatively few 
anti-protons in the dataset. As a result, we use a less re¬ 
strictive PID cut for the anti-proton, keeping events with 
Atof between -1.8 ns and 1.0 ns. Because of the absence 
of background reactions that produce p, this PID cut is 
the most stringent requirement in this channel. We op¬ 
timize cuts on the and MM{K^) observables in 

the same method as for the A —)■ ml channels, yielding 
optimal cut widths of wi = 5.00 x 10“^ GeV^/c"*^ and 
W 2 = 4.25 X 10“^ GeV/c^, respectively, and a signal effi¬ 
ciency of e = 4.98%. 


C. A ^ K%v 


In addition to the charged decay modes, we search for 
the decay of A —)• K'gV, using the dominant charged decay 
mode ATg —>■ tt+tt”. (Observation of the —>■ 27r selects 
ATg rather than Ar£.) Because of this reaction’s final- 
state neutrino, we do not have access to /W(A'g,z/); 
thus, the analysis described for the charged decay modes 
is not appropriate for this channel. In addition, the un¬ 
measured momentum of the v limits our analysis con¬ 
straints and we can expect more background to pass the 
optimized cuts. 

We begin the background separation process by apply¬ 
ing two-dimensional PID cuts based on Atof and mo¬ 
mentum to the charged final state particles (recoil Ar+ 
and the from A'° decay) similar to those for the 
A —7r+/i“ channel (see Fig. [^. We then optimize a 
two-dimensional cut motivated by the particulars of this 
decay. The first is a symmetric cut on centered 

at 0 (the mass of the v is negligible) with width wi. 

The second cut identifies pairs that are produced 
from Ar° decay by inspecting the tt* opening angle, {0'^), 
in the c.m. frame: 


= cos 


-1 


( P+-P- \ 

ViT+iiT-iy ’ 


( 11 ) 


where p± are the momenta of the decay pions indexed by 
charge. Due to the break-up energy associated with the 
—>■ TT+TT decay, 0^ is constrained to a narrow range 
for a given value of momentum. However, for 7r“'"7r“ 
pairs that do not come from decay, we expect only 
the constraints associated with momentum conservation 
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FIG. 4. (Color Online) Atof vs p for A —> fi~ Monte Carlo. The red curves define the boundaries of the PID cuts for each 
particle type. Shown are the recoil (left), decay 7r+ (center), and decay /r“ (right). The bands with Atof > 1.0 ns in the 
kaon plot represent from kaon decays. 


for the entire event, i.e. much less correlation between 
the pion momenta. Distributions of 9'^ vs magnitude of 
momentum, px, for data and signal MC are shown 
in FiglU We separate —>■ tt+tt events by making a 

two-dimensional cut on 9^ vs px- To obtain a description 
of the correlation between 0^ and px, we fit the two- 
dimensional histogram of the two observables and found 
adequate description with the function 


0 459 

f{px)= 0.319- L + 


2.23 


0.641 


Px+0.5 ' (px+O.^)"^ {px+0.5)^' 

( 12 ) 

The cut width, W 2 , is implemented by keeping events for 
which 


{1-W2)f{px)<9^<{l+W2)f{px). (13) 


The optimization process tests 1600 pairings of cut 
widths with wi S [0,0.15] GeV^/c^ and u ;2 G [0,0.2]. Be¬ 
cause / NV{K^, v) is not accessible, we use the MM ) 
distributions to estimate €{C) and B{C), and ultimately 
to identify the signal. We choose the blinded signal 
region to be within 0.03 GeV/c^ of the A mass peak 
in the MM{K'^) spectrum determined from standard 
model data and MC, 1.186 GeV/c^, and the side-band 
regions to be within 0.15 GeV/c^ of the peak (exclud¬ 
ing the signal region). We find the optimal widths to be 
wi = 1.875 X 10“^ GeV^/c"^ and W 2 = 0.0600, yielding 
an efficiency of 2.23% and 239.25 estimated background 
events in the signal region. The MM{K^) distribution 
for data events after unblinding is shown in Fig. 


D. Assessing background signatures 


Because of the sensitive nature of a positive signal 
identification, understanding the signature (shape) of the 
background in signal identification histograms is cru¬ 
cial. If a peak is present when the data in unblinded, 
we must be sure that the peak represents BNV signal, 
and is not merely an unfortunate distortion of the back¬ 
ground events due to the cuts used. In order to assess 
the signature of the background events in INV{A,B) 
and MM{K^) distributions, we generated 5 x 10^ Monte 
Carlo events for each of six background reactions: 


• 7P —>■ K'^A —>■ K~^pTT~ 

• "fp ^ K^TP —>■ K~^p'n~^ 

• 7P —>■ pe^e~ 

• 7P —>■ ppA pL~ 

• JP ^ PTT^'K~ 

• 7 P —>■ pK^ K~ 

• 7P —>■ p7r+7r“7r*^ 


These reactions were chosen either for their abundance 
in the dataset (combination of large cross section and 
detectability in CLAS) or for their similarity to the BNV 
channels investigated (similar final-state particles). After 
applying the optimized cuts for each BNV channel to 
these MC events, we found that very few events pass the 
cuts; no channel’s cuts allow more than one background 
MC event into the signal region of the INV{A,B) or 
MM{K^) histograms. Thus, we claim that none of the 
background reactions investigated create an excess in the 
signal regions. 
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FIG. 5. (Color Online) Plots of P vs cut widths wi (GeV^/c'*) and W 2 (GeV/c^) (a), data INV{tt^ fi~) (b), Monte Carlo 
MM{K'^) (c), and Monte Carlo MM^ (d) for the A —>■ 7r+p“ channel. In (b), the vertical red lines show the boundaries of the 
blinded signal region, and the green vertical lines show the boundaries of the side-band regions. In (c) and (d), the grey and 
red lines show the center and boundaries of the optimal cuts, respectively. 


V. RESULTS to the unblinded data. For the nine decay modes where 

the final state can be completely reconstructed, we found 
A. Charged decays the number of observed events in the signal region, Noh&, 

to be between 0 and 2, consistent with background esti- 

After the selection criteria are finalized using the 
Monte Carlo and side-band studies, we applied these cuts 
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p (GeV/c) 


FIG. 6. (Color Online) Kg —>■ tt+tt opening angle vs p for A 
red curves demonstrate the boundaries of the optimized cut. 



FIG. 7. (Color Online) Unblinded MM(K+) distribution for 
data events after application of cuts optimized for the A —>■ 
KgV channel. The red vertical lines show the boundaries of 
the blinded signal region, and the green vertical lines show 
the boundaries of the side-band regions. 


mates from the cut optimization studies. For these decay 
modes we used the Feldman-Cousins approach 121! to de¬ 
termine upper limits on the reconstructed signal yields, 
iVuL. 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 


p (GeV/c) 

Kgv Monte Carlo (left) and data (right) after PID cuts. The 


The Feldman-Cousins approach provides a way to esti¬ 
mate upper confidence limits for null results. The inputs 
are the expected number of background events (A"eb) and 
the observed number of events (A^obs)- We estimated the 
expected number of background events from the side¬ 
bands in the data (see previous section). iVobs and TVeb 
values for each decay mode are shown in Table |TTj Frac¬ 
tional numbers are given when only one event was ob¬ 
served in a side-band region that spans a greater range 
than the signal region. 

These provided the input to the Feldman-Cousins 
method and we quote the upper limit on the recon¬ 
structed signal yield at 90% confidence level (A^ul), 
shown in Table |TT] To calculate the upper limit on the 
branching fraction (Sul), we used the efficiency (e) for 
each decay mode as determined from Monte Carlo stud¬ 
ies and the total number of jp AT+A —J, 
events produced during the data-taking period, Afp^od 
(see Eqn. |^: 


(14) 

cTy prod 


where BpT^ = 0.639 ± 0.005 is the A ^ pn branching 
fraction m- The Sul values for all charged decay modes 
are shown in Table IE 
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TABLE 11. Inputs to and final calculations for the upper limits on the branching fractions. wi and W 2 are the optimized 
cut widths for each channel. For the charged decay modes, Wi and W 2 are the widths of the cuts on MM^ (in GeV^/c"^) 
and MM{K^) (in GeV/c^), respectively. For A —>■ KgU, wi and W 2 are the widths of the cuts on (in GeY^/c^) and 

Kg —^ 7r"'"7r“ opening angle, respectively, e is the efficiency for each reaction, calculated using Monte Carlo studies described 
in the text, Aeb is the number of expected background events, calculated from side-band studies and confirmed with Monte 
Carlo studies, Aobs is the number of observed events in the data, A^ul is an upper limit on the number of signal events in the 
data, calculated using the Feldman-Cousins technique for the charged decay modes and likelihood scanning technique for 
the A —> KgU channel, and Bul is the upper limit on the A branching fraction for each decay mode. 


Reaction 

Wi 

W2 

£(%) 

Aeb 

Aoba 

Aul 


Bul 

A ^ K+e- 

2.50 X 10“^ 

0.01625 

4.13 

0 

1 

4.36 

2 

X 10"“ 

A-> A+p- 

3.25 X 10"^ 

0.0125 

4.42 

0 

2 

5.91 

3 

X 10"® 

A^ A-e+ 

1.80 X 10"® 

0.01375 

4.63 

0 

1 

4.36 

2 

X 10"® 

A-> 

3.00 X lO"'^ 

0.0300 

4.40 

0 

2 

5.91 

3 

X 10"® 

A —>■ 7r"'"e“ 

2.75 X 10""^ 

0.00900 

7.02 

0 

0 

2.44 

6 

X 10"'^ 

A — >■ pT 

3.25 X lO"'^ 

0.00900 

7.91 

0 

0 

2.44 

6 

b- 

1 

O 

r—H 

X 

A —>■ 7r“e^ 

4.75 X lO""^ 

0.0125 

8.65 

0.75 

0 

1.94 

4 

X 10"^ 

A —> TT~ 

3.50 X lO"'^ 

0.00900 

7.92 

0.25 

0 

2.44 

6 

b- 

1 

o 

r—H 

X 

A —>■ pvT'*’ 

5.00 X lO""^ 

0.0425 

4.98 

0 

0 

2.44 

9 

X 10"'^ 


0.01875 0.0600 2.23 239.25 -3.88 14.1 2 x 10"® 


B. A ^ K°sv 

For the K^v decay mode, a BNV signal would manifest 
itself as a peak in the MM{K'^) distribution at the A 
mass. When we unblinded the MM{K'^) histogram, we 
observed no such peak and found a number of events in 
the signal region that is consistent with the background 
study above. The number of events in the signal region 
is much larger than is normally handled by the Feldman- 
Cousins approach; thus, we perform a likelihood scan to 
determine the upper limit on Aul- 

We performed an unbinned maximum likelihood fit to 
the data in this region using an exponential function 
probability density function (PDF) to describe the back¬ 
ground and a sum of two Gaussians to describe the sig¬ 
nal. The shape of the background was allowed to vary 
in the fit, as are the numbers of signal and background 
events. The parameters describing the Gaussians (be. 
means and widths) were fixed to values determined from 
Monte Carlo studies. 

The fit converged to a central value of Aobs = — 3.88± 
8.9 signal events, consistent with 0 signal events. To 
check whether or not this negative value is of con¬ 
cern, we sampled from a distribution described by the 
background parameters returned by the fit to generate 
1000 mock “background-only” samples and fit them to a 
background-plus-signal hypothesis. About 50% of these 
fits returned a negative value for the signal and about 
35% returned a value more negative than what was found 
in the data. We determined that the negative value is an 
artifact of fitting to a small number of points using a 
function with as much freedom as we use. We note that 
nowhere does the total PDF go negative. 

To calculate an upper limit on the signal yield, we 
scaned the likelihood function by performing a series of 
fits where the signal yield (A) is varied around the best fit 
value Aobs and the other parameters were refit to map out 


the difference in the In-likelihood: Aln£ = In£(Aobs) ~ 
ln£(A). We integrated the function y = over 

A. We ignored the unphysical region with A < 0 and 
calculate the integral for A > 0. We note the value of A 
which encloses 90% of the area above A = 0 and interpret 
this as the upper limit on the signal yield returned by 
the fit at 90% confidence level. This procedure returns 
an upper limit (Ay^) of 14.1 signal events. 


C. Experimental uncertainties 

Uncertainty in ^Bul comes from the world average of 
the A —>■ p7r+ branching fraction (0.8%) [T^] and sta¬ 
tistical and systematic uncertainties from the extraction 
of Aprod and cut efficiencies for each BNV channel. We 
found a 6.1% relative uncertainty in Aprod by combining 
the 0.02% systematic uncertainty due to the INV{p, tt~) 
peak fitting procedure, the « 6% systematic uncertainty 
in GLAS acceptance calculation (taken from previous hy- 
peron production analysis m), and the 0.64% statistical 
uncertainty in (estimated with binomial statistics). 
We estimated the uncertainty in e for each BNV chan¬ 
nel by comparing the effects of optimized MM{K^) and 
cuts on MG and standard-model data distribu¬ 
tions, and found it to be « 7.6%. We combined all un¬ 
certainties in quadrature to find a relative uncertainty in 
Bul of « 9.8%. With this estimate of the combined un¬ 
certainty in hand, we quote the final Bul results to one 
significant figure (see Table [lT| . 


VI. SUMMARY 

The analysis described here represents the first search 
for baryon- and lepton-number violating decays of the A 
hyperon. Though similar studies have been performed 
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with much higher sensitivities for decays of the nucleon, 
this study offers the first direct probe of BNV processes 
involving strange quarks in the initial state. Using a 
dataset for photoproduction off of the proton collected 
with the CLAS detector at Jefferson Laboratory con¬ 
taining roughly 1.8 x 10® reconstructed Standard Model 
A —> pTT~ decays, we have searched via blinded analysis 
for BNV decays of the A to either meson-lepton pairs or 
to p7r+. We found no BNV signal in any of the ten decay 
channels investigated, and set upper limits on branching 
fraction for each of the processes studied in the range 
7 X lO-"^ to 2 X 10"®. 
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